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In these lectures we provide an introduction to the theory of QCD at very high baryon density. 
We begin with a review of some aspects of quantum many-body system that are relevant in the QCD 
context. We also provide a brief review of QCD and its symmetries. The main part of these lectures 
is devoted to the phenomenon of color superconductivity. We discuss the use of weak coupling 
methods and study the phase structure as a function of the number of flavors and their masses. We 
also introduce effective theories that describe low energy excitations at high baryon density. Finally, 
we use effective field theory methods in order to study the effects of a non-zero strange quark mass. 
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I. INTRODUCTION 



In these lectures we wish to provide an introduction to recent work on the phase structure of QCD at non-zero 
baryon density. This work is part of a larger effort to understand the behavior of matter under "extreme" conditions 
such as very high temperature or very large baryon density. There are several motivations for studying extreme QCD: 

• Extreme conditions exist in the universe: About 10~^ sec after the big bang the universe passed through a state 
in which the temperature was comparable to the QCD scale. Much later, matter condensed into stars. Some of 
these stars, having exhausted their nuclear fuel, collapse into compact objects called neutron stars. The density 
at the center of a neutron star is not known very precisely, but almost certainly greater or equal to the density 
where quark degrees of freedom become important. 

• Exploring the entire phase diagram is important to understanding the phase that we happen to live in: We cannot 
properly understand the structure of hadrons and their interactions without understanding the underlying QCD 
vacuum state. And we cannot understand the vacuum state without understanding how it can be modified. 

• QCD simplifies in extreme environments: At scales relevant to hadrons QCD is strongly coupled and we have 
to rely on numerical simulations in order to test predictions of QCD. In the case of large temperature or large 
baryon density there is a large external scale in the problem. Asymptotic freedom implies that the bulk of the 
system is governed by weak coupling. As a result, we can study QCD matter in a regime where quarks and 
gluons are indeed the correct degrees of freedom. 

• Finally, extreme QCD tries to answer one of the simplest and most straightforward questions about the behavior 
of matter: What happens if we take a piece of material and heat it up to higher and higher temperature, or 
compress to larger and larger density? 

There are several excellent text books and reviews articles that provide an introduction to QCD and hadronic 
matter at finite temperature Q, |^ Q . In these lectures we will focus on matter at high baryon density and small 
or zero temperature. This is the regime of the "condensed matter physics" of QCD Q. Ordinary condensed matter 
physics is concerned with the overwhelmingly varied appearance and rich phase diagram of matter composed of 
electrons and ions. All phases of condensed matter ultimately derive their properties from the simple laws of quantum 
electrodynamics. We expect, therefore, that the simple laws of QCD will lead to a phase diagram of comparable 
diversity. In fact, since there is only one kind of electron, but several flavors and colors of quarks, we might expect 
new and unusual phases of matter never before encountered. 

These lectures are organized as follows. In sections IIIIIVI we review a number of simple many body systems that are 
relevant to the behavior of QCD matter in different regimes. In order to keep the presentation simple, and to make 
contact with well-known properties of other many body system, we phrase our discussion not in terms of quarks and 
gluons, but in terms of generic fermions and bosons interacting via short range range forces. In section^we provide 
a brief introduction to QCD and its symmetries. Sections IVIIIIXI form the main part of these lectures. We introduce 
the phenomenon of color superconductivity, study the phase structure in weak coupling, and introduce effective field 
theories that allow systematic calculations of the properties of dense QCD matter. Other aspects of high density 
QCD are discussed in the many excellent reviews on the subject 0,111113. 



II. FERMI LIQUIDS 
A. Introduction 

In this section we wish to study a system of non-relativistic fermions interacting via a short-range interaction d |^ . 
The lagrangian is 

Co^yJ^(^^o + ^^i^-^{i:HT■ (1) 

The coupling constant Cq is related to the scattering length, Cq — Aiva/m. Note that Co > corresponds to a 

repulsive interaction, and Co < is an attractive interaction. The lagrangian equ. is invariant under the U{1) 
transformation i/; — > e"^-)/). The U{V) symmetry implies that the fermion nmnbcr 

N ^ ( d^xV'V (2) 
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b) 






FIG. 1: Leading order Feynman diagrams for the ground state energy of a dilute gas of fermions interacting via a short range 
potential. 



is conserved. As a consequence, it is meaningful to study a system of fermions at finite density p — N/V. We will do 
this in the grand-canonical formalism. Wc introduce a chemical potential /i conjugate to the fermion number N and 
study the partition function 



-l3{H-iiN) 



(3) 



Here, H is the Hamiltonian associated with C and /3 = l/T is the inverse temperature. The trace in equ. Q runs over 
all possible states of the system, including all sectors of the theory with different particle number A^. The average 
number of particles for a given chemical potential jj. and temperature T is given by {N) = T {d log Z) / {d p.) . At zero 
temperature the chemical potential is the energy required to add one particle to the system. 

There is a formal resemblance between the partition function equ. Q and the quantum mechanical time evolution 
operator U = cxp{—iHt). In order to write the partition function as a time evolution operator we have to identify 
(3 it and add the term —jiN to the Hamiltonian. Using standard techniques we can write the time evolution 
operators as a path integral 0] 



Z 



Here, Ce Is the euclidean lagrangian 



D^Dip^ exp 



IpUdr - p- 
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dr / d xCe 



2m 



(4) 



(5) 



The fermion fields satisfy anti-periodic boundary conditions ip{l3) = —'ip{0). Equation Q is the starting point of the 
imaginary time formalism in thermal field theory. The chemical potential simply results in an extra term —fi^'^ip in 
the lagrangian. From equ. lO we can easily read off the free fermion propagator 



5a 



2m 



(6) 



where a, (3 are spin labels. We observe that the chemical potential simply shifts the four-component of the momentum. 
This implies that we have to carefully analyze the boundary conditions in the path integral in order to fix the pole 
prescription. The correct Minkowski space propagator is 



Oal3 



pa 



+ iSsgn{ep] 



= Sa 



Q{p-Pf) , Q{pF-p) 



Po 



iS 



Po 



iS 



(7) 



where ep = Ep — /i, Ep = p"^ /(2m) and 5 — s- 0^. The quantity pp = •\/2m/^ is called the Fermi momentum. We 
will refer to the surface defined by the condition \p\ — pf as the Fermi surface. The two terms in equ. ((TJ have 
a simple physical interpretation. At finite density and zero temperature all states with momenta below the Fermi 
momentum are occupied, while all states above the Fermi momentum are empty. The possible excitation of the system 
are particles above the Fermi surface or holes below the Fermi surface, corresponding to the first and second term in 
equ. (|7|). The particle density is given by 



P 



d^p 



0+ 



d^p 
(2^ 



Q{PF 



(8) 



FIG. 2: Leading order Feynman diagrams that contribute to the photon polarization function in a non-relativistic Fermi liquid. 
The tadpole diagram shown in the right panel only appears in the spatial part of the polarization tensor. 



As a first simple application we can compute the energy density as a function of the fermion density. For free fermions, 
we find 



d^P p., ^ ^ Pf 

^E,e{pF~p)^-p—. 



I 7^E,eipF~p)^-p^. (9) 



We can also compute the corrections to the ground state energy due to the interaction ^Co{ip^tp)'^ . The first term is 
a two- loop diagram with one insertion of Cq, see Fig. ^ We have 



r,3 
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We should note that equ. H10() contains two possible contractions, usually called the direct and the exchange term. If 
the fermions have spin s and degeneracy g = (2s + 1) then equ. has to be multiphed by a factor g{g — l)/2. We 
also note that the sum of the first two terms in the energy density can be written as 

£-P^(l + liPFa) + ..), (11) 



2m \5 3tt 

which shows that the Co term is the first term in an expansion in ppa, suitable for a dilute, weakly interacting, Fermi 
gas. The expansion in {pfo) was carried out to order {pFa)"^ by Huang, Lee and Yang 0, Since then, the 
accuracy was pushed to ©((pi^a)"* log(p_Fa)) see for a modern perspective. The effective lagrangian can also 
be used to study many other properties of the system, such as corrections to the fermion propagator. Near the Fermi 
surface the propagator can be written as 

"'^ Pa-VFi\p\-pF)+iSsgn{\p\-pF)' 

where Z is the wave function renormahzation and vf — PF/fn* is the Fermi velocity. Z and m* can be worked out 
order by order in {pfo), see The main observation is that the structure of the propagator is unchanged even 

if interactions are taken into account. The low energy excitations are quasi-particles and holes, and near the Fermi 
surface the lifetime of a quasi-particle is infinite. This is the basis of Fermi liquid theory, see Sec. |nj We should note, 
however, that for nuclear systems the (pfo-) expansion is not particularly useful since the nucleon-nucleon scattering 
length is very large. Equ. pi(l is of interest for trapped dilute Fermi gases. 



B. Screening and damping 

An important aspect of the dilute Fermi gas is the response to an external electromagnetic field. As a simple 
example we will consider the case of an a static electric field. The coupling of the gauge field is given by e^oV'^V'- 
The medium correction to the photon propagator is determined by the polarization function 

noo(9) = y d4a;e"^«-(V'V(0)V'V(x)). (13) 

The one-loop contribution is given by 

n ( ) - 2 / '^^^ ^ ^ (14) 

°° J (27r)4 q^+pQ- ep+q + i(5sgn(ep+,) Po - f-p + i&gn(ep) ' 
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Performing the po integration by picking up the pole we find 



where we have introduced the Fermi distribution function Up = Q{pf ~ p)- We observe that in the Umit (f 
the polarization function only receives contributions from particle-hole pairs that are closer and closer to the Fermi 
surface. On the other hand, the energy denominator diverges in this limit because the photon can excite particle-hole 
pairs with arbitrarily small energy. As a result we get a finite contribution 

noo(<Zo = 0,,-0)^e y (^5^=e^, (16) 

which is proportional to the density of states on the Fermi surface. Equ. H16|l implies that the static photon propagator 
in the limit (f ^ is modified according to l/ff^ l/('f^ + where 

2 2fPFm\ 

is called the Debye mass. The factor N — (ppm)/ {2tt'^) is equal to the density of states on the Fermi surface. In a 
relativistic theory we find the same result as in equ. H17() with the density of states replaced by the correct relativistic 
expression N ~ {ppEp)/ {2tt'^). The Coulomb potential is modified as 

V{r) = -e , (18) 

r 

where ro = ^/niD is called the Debye screening length. The physics of screening is very easy to understand. A test 
charge can polarize virtual particle-hole pairs that act to shield the charge. 

In the same fashion we can study the response to an external vector potential A. The coupling of a non-relativistic 
fermion to the vector potential is determined in the usual way by replacing p ^ p + eA. Since the kinetic energy 
operator is quadratic in the momentum wc find a linear and a quadratic coupling of the vector potential. The one-loop 
diagrams that contribute to the polarization tensor are shown in Fig. [5] In the limit of small external momenta we 
find 



n,, (g) = e^ml / ^ v,v, ""^'^ ^'^^^ - , (19) 



dQ ( vkiq ■ p) \ y 

Att { 'qo-vkiq-p) 3 

where v — p/m is the Fermi velocity. In the limit qo = the polarization tensor vanishes. There is no screening of 
static magnetic fields. For non-zero qo the trace of the polarization tensor is given by 

n^i{q) ^mjj—logl ■ . (20) 



2q \qo + vq 

The result equ. (|20|l has an imaginary part if vq > qo. This phenomenon is known as Landau damping. The photon 
is loosing energy to electrons in the Fermi liquid. For a discussion in the context of kinetic theory we refer the reader 
to [3. 



III. BOSE CONDENSATION 



In this section we introduce some general features of bosonic systems at finite density. We will consider a charged 
relativistic boson described by the Lagrange density 

Co = id^rM.cl>) ~ rn^rcb - X{<f>*<jyf. (21) 

Note that A has to be positive in order for the theory to be well defined. This corresponds to a repulsive interaction 
between the bosons. The lagrangian has a U{\) symmetry (j) — > e~''^(j3. The corresponding conserved charge is 

Q [ d^xi{(t)*do(j)- (t)dQ(l)*) . (22) 
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FIG. 3: Spectrum and charge density for a charged scalar field as a function of the chemical potential. 



Note that the charge density p contains not only the field but also the canonically conjugate momentum d^cj). 
This means that the chemical potential modifies the integration over the canonical momenta in the path integral 
representation of the partition function. The resulting Minkowski space path integral is 



D(f)D(f)* exp J (fxCj , (23) 
with 

C = {do+ iAi)0*(9o - - (V0*)(V0) - m^^*^ ~ K<P*<pf- (24) 

There is a simple argument that fixes the form of the lagrangian equ. I|24|l. The argument is based on the observation 
that we can promote the global U{1) symmetry to a local symmetry by adding a f/(l) gauge field to the lagrangian. 
The charge density is obtained by varying the effective action with respect to the gauge potential. This implies that 
the chemical potential has to enter the lagrangian like the time component of a gauge field. 

We can study the effect of a chemical potential in the mean field approximation. The classical effective potential 
for the field (j) is given by 

V{4>) = {m^-p^Wcl^) + \{<ly*<f)\ (25) 

For fi > m the quadratic term is positive and the minimum of the effective potential is at (0) = 0. For /i > m the 
origin is unstable and 

7 2 

{4'? = (26) 

This state is a Bose condensate. The charge density is 

p=^(M^-m2). (27) 

In a non-interacting Bose gas the chemical potential cannot be larger than the mass of the boson. In our case, repulsive 
interactions limit the growth of the density and the chemical potential can take any value. We can also compute the 
spectrum as a function of the chemical potential. We write </> = {(j)) + Xi + *X2 and expand the effective action to 
second order in xi,2- For fi < m we find two modes with energies E{p = 0) — m± fi. Bose condensation sets in when 
the lower mode reaches zero energy. Above the onset of Bose condensation we find 

Ei{p^O) = 0, E2{p = Q) = V^f^^ - 2^2. (28) 

Bose condensation breaks the U{1) symmetry spontaneously and the spectrum contains one Goldstone boson. It is 
also interesting to study the dispersion relation of the Goldstone mode in more detail. For small momenta we find 
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ZS ZS' BCS 

FIG. 4: Second order diagrams that contribute to particle-particle scattering. The three diagrams are known as ZS (zero 
sound), ZS' and BCS (Bardeen-Cooper-Schrieffer) contribution. 



This shows that at the phase transition point the velocity of the Goldstone mode is zero. Far away from the transition 
the velocity approaches v = c/VS. Bose condensates have many remarkable properties, most notably the fact that 
they can flow without viscosity. These properties can be derived from the effective action for the Goldstone mode. It 
was shown, in particular, that this effective action is equivalent to superfluid hydrodynamics [TsLIT^. 



IV. SUPERCONDUCTIVITY 



A. BCS instability 



One of the most remarkable phenomena that take place in many body systems is superconductivity. Superconduc- 
tivity is related to an instability of the Fermi surface in the presence of attractive interactions between fermions. Let 
us consider fermion-fermion scattering in the simple model introduced in Sect. ^] At leading order the scattering 
amplitude is given by 

^al3js{Pl,P2,P3,P4) = Co {Saj6f3S - SasSfjj) ■ (30) 

At next-to-leading order we find the corrections shown in Fig. 0] A detailed discussion of the role of these corrections 
can be found in |1, The BCS diagram is special, because in the case of a spherical Fermi surface it can lead 
to an instability in weak coupling. The main point is that if the incoming momenta satisfy pi ~ —p2 then there are 
no kinematic restrictions on the loop momenta. As a consequence, all back-to-back pairs can mix and there is an 
instability even in weak coupling. 

For pi — —p2 and Ei = E2 = E the BCS diagram is given by 

^ap-iS ^ Cl {5a-^5i3S - SasSp-i) ( , '^4 TT-; "^7 ' (^1) 

J (27r)^ E + qa- eq + idsgn(e,) E - qo ~ eq + zdsgn(e,) 

The loop integral has an infrared divergence near the Fermi surface as -E — ^- 0. The scattering amplitude is proportional 
to 

where Eq is an ultraviolet cutoff. Ea u. 13 211 can be interpreted as an effective energy dependent coupling that satisfies 
the renormalization group equation [13, UM 



j^dCo _ 2 



with the solution 



^"("^^ - l + NCo?E!)tgiEo/Ey 
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where A'' = {ppm) / {2tt'^) is the density of states. Equ. (|34|l shows that there are two possible scenarios. If the initial 
coupling is repulsive, Cq{Eq) > 0, then the renormalization group evolution will drive the effective coupling to zero 
and the Fermi liquid is stable. If, on the other hand, the initial couphng is attractive, Cq{Eq) < 0, then the effective 
coupling grows and reaches a Landau pole at 

Ecrit ~ Eq exp ( -^rTTTT-TTTTT ) • (35) 



N\CoiEo)\ 

At the Landau pole the Fermi liquid description has to break down. The renormalization group equation does not 
determine what happens at this point, but it seems natural to assume that the strong attractive interaction will lead 
to the formation of a fermion pair condensate. The fermion condensate {e^^ipaipp) signals the breakdown of the U{1) 
symmetry and leads to a gap A in the single particle spectrum. 

The scale of the gap is determined by the position of the Landau pole, A ^ Ecrit- A more quantitative estimate of the 
gap can be obtained in the mean field approximation. In the path integral formulation the mean field approximation 
is most easily introduced using the Hubbard-Stratonovich trick. For this purpose we first rewrite the four-fermion 
interaction as 

^(^t^)2 ^ ^ {{i;^a2tP^)(tPa2tP) + (V'V2aV^)(^Aaa2V)} , (36) 

where we have used the Fierz identity 2S°''^S^p = S^pS""^ + {a)°'P{a)'<f^ . Note that the second term in equ. (|36|l vanishes 
because ((72^) is a symmetric matrix. We now introduce a factor of unity into the path integral 

[ DAexpf^] , (37) 



Za J \ Co 

where we assume that Co < 0. We can eliminate the four-fermion term in the lagrangian by a shift in the integration 

variable A. The action is now quadratic in the fermion fields, but it involves a Majorana mass term ipa2Atfj + h.c. The 

Majorana mass terms can be handled using the Nambu-Gorkov method. We introduce the bispinor = {ijj , (T2) 
and write the fermionic action as 

S=l f-^^f(Po-/p ^ U. (38) 

2 J (27r)4 \ A* po + epj ' 



Since the fermion action is quadratic we can integrate the fermion out and obtain the effective lagrangian 

lTr[log(Go-iG)]+^ 



/: = lTr[log(Go-iG)]+-^|Ap, (39) 



where G is the fermion propagator 



The diagonal and off-diagonal components of G{p) are sometimes referred to as normal and anomalous propagators. 
Note that we have not yet made any approximation. We have converted the fermionic path integral to a bosonic one, 
albeit with a very non-local action. The mean field approximation corresponds to evaluating the bosonic path integral 
using the saddle point method. Physically, this approximation means that the order parameter does not fluctuate. 
Formally, the mean field approximation can be justified in the large N limit, where N is the number of fermion fields. 
The saddle point equation for A gives the gap equation 



Performing the po integration we find 



1^^ fj^ 

2 J (2^)3 ^fT^2 



1-^/7^^=^- (42) 



Since ep = Ep — fj, the integral in equ. (|42|l has an infrared divergence on the Fermi surface \p\ ^ pp- As a result, 
the gap equation has a non-trivial solution even if the coupling is arbitrarily small. The magnitude of the gap is 
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A ~ Aexp(— 1/(1 Col iV)) where A is a cutoff that regularizes the integral in equ. H42() in the ultraviolet. If we treat 
equ. (P) as a low energy effective field theory we should be able to eliminate the unphysical dependence of the gap on 
the ultraviolet cutoff, and express the gap in terms of a physical observable. At low density, this can be achieved by 
observing that the gap equation has the same UV behavior as the Lipmann-Schwinger equation that determines the 
scattering length at zero density 

I^^/t^I^- (43) 



A-Ka 2 J {2ttY Ep 

Combining equs. (|42|) and H43|) we can derive an UV finite gap equation that depends only on the scattering length, 

I f d^p r 1 



47ra 2 J (27r)3 L ^ /^J^^ 

A careful analysis gives [T^l20j 



V 2pF\a\J 

For neutron matter the scattering length is large, a — —18.8 fm, and equ. (|45|l is not very useful, except at very small 
density. Calculations based on potential models give gaps on the order of 2 MeV at nuclear matter density. 

In the limit of very high density we can eliminate the cutoff dependence using a method introduced by Weinberg [2lj . 
Weinberg defines a renormalized effective potential and shows that the renormalization scale dependence of the effective 
potential is canceled by the scale dependence of the coupling. The effective coupling satisfies the renormalization group 
equ. The gap is determined by the effective coupling at the energy scale Eq. In practice, this would typically 

be the energy scale at which the four-fermion interaction is matched against a more microscopic description in terms 
of meson (nuclear physics) or phonon exchange (condensed matter physics) . 



B. Fermi liquid, revisited 

Our discussion of Fermi liquids in Sect. IITI and in the previous section was based on the simple model defined in 
equ. JQ). In this section we shall briefly discuss the structure of fcrmionic many-body systems in the case of more 
general interactions. We will restrict ourselves to systems that can be described in terms of purely fcrmionic actions, 
with all other degrees of freedom integrated out. For more details we refer the reader to | il7j . 

We can view the model defined by equ. 0]) as an example of an effective field theory, valid for momenta close to 
the Fermi surface. In order to construct an effective field theory we have to write all possible interactions that are 
allowed by the symmetries of the theory. The effective action of rotationally invariant, non-relativistic Fermi system 
is given by 

i^HP4:)'4'HP3)ll'{P2)ll'{Pl)U{p4,P3,P2,Pl), (46) 

where vp = dtp/ {dp) and Ip — \p\ —pp. We have suppressed the spin indices of the potential U . The power counting for 
the effective theory can be established by studying the scaling behavior of all allowed operators under transformations 
of the type Ip — > sip that scale the momenta towards the Fermi surface. Writing tp = vplp + 0{lp) we see that as 
s — > only the Fermi velocity survives, the detailed form of the dispersion relation is irrelevant. Using this method 
we can also see that with the exception of special kinematic regimes the four-fermion interaction is irrelevant. We 
already saw that one exception is provided by the BCS interaction 

U{-p3,P3,-Pi,Pi) = V{pi -ps) = ^ViPiipi ■P3), (47) 

I 

where Pi{x) are Legendre polynomials. At tree level V{x) is a marginal operator, that means it is invariant under 
rescaling the momenta towards the Fermi surface. This changes at one-loop level. If any of the couplings Vi is attractive 
then this coupling will grow according to the renormalization group equ. H33|l and eventually reach a Landau pole. 
If there is more than one attractive coupling Vi then the ground state is determined by which coupling reaches its 
Landau pole first. If all Vi are repulsive then the BCS potential becomes irrelevant as the evolution approaches the 
Fermi surface. 



S = 



n 
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BCS Landau 

FIG. 5: Kinematic configurations for fermion-fermion scattering that correspond to four-fermion operators that are marginal 
at tree level. The left panel shows BCS (back-to-back) scattering and the right panel shows forward scattering. 



In this case there is another kinematic regime that becomes important. We can take any two momenta on the Fermi 
surface, not necessarily back-to-back, and find an allowed final state. Energy and momentum conservation implies 
that pi ■ P2 = P3 ■ Pa- In two-dimensions this would restrict the scattering to be either forward or exchange, but in 
three dimensions there is a circle of allowed final states parametrized by the angle 012,34 between the planes spanned 
by the incoming and outgoing momenta, see Fig. [S] The interaction is 

UiP4,P3,P2,Pl)\p^.p^=p^.p^ =F{pi -^2, 012,34) (48) 

The function F{x, 0) is called the Landau function and its Legendre coefficients are referred to as Landau parameters. 
The Landau parameters remain marginal at one-loop order. A many body system characterized by vp and Fi is called 
a Landau Fermi liquid |2^|2^. A Landau liquid behaves qualitatively like a dilute, weakly interacting, Fermi liquid, 
even if the interaction is not weak and the system is not dilute. In particular, the excitations of a Landau liquid are 
quasi-particles and holes. The Landau parameters encode the quasi-particle interaction and can be used to compute 
observables like the equation of state and the response to external fields. 

One can show that operators involving more than four fermion fields are irrelevant near the Fermi surface. This 
does not imply that n > 4 fermion operators play no role at all. For example, if one of the Vi has a Landau pole at 
energy Eq then the six fermion interaction still has a finite coupling at this scale and will cause observable eflfects, see 
Sec. IIX El for an example. Also, just because the Vi are the only operators that cause instabilities in weak coupling 
does not imply that other operators cannot have instabilities in strong coupling. For example, dilute nuclear matter 
may have a phase characterized by alpha particle condensation rather than superconductivity. 



C. Landau-Ginzburg theory 



In this section we shall study the properties of a superconductor in more detail. For definiteness, we will consider 
a system of electrons coupled to a U{1) gauge field A^^. The order parameter <& — (e^^ipaipfj) breaks U{1) invariance. 
Consider a gauge transformation 

A^^A^ + d^A. (49) 

The order parameter transforms as 

$ exp(2ieA)$. (50) 

The breaking of gauge invariance is responsible for most of the unusual properties of superconductors 119 • This 
can be seen by constructing the low energy effective action of a superconductor. For this purpose we write the order 
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parameter in terms of its modulus and phase 



$(a:) = cxp(2ie(/)(x))$(x). 

The field (f) corresponds to the Goldstone mode. Under a gauge transformation (l){x) 
restricts the form of the effective Lagrange function as 



L = - 



(51) 

(x) + A{x). Gauge invariance 



(52) 



There is a large amount of information we can extract even without knowing the explicit form of Lg. Stability implies 
that Afj^ ~ dfj_(f) corresponds to a minimum of the energy. This means that up to boundary effects the gauge potential 
is a total divergence and that the magnetic field has to vanish. This phenomenon is known as the Meissner effect. 

Equ. (|52|l also implies that a superconductor has zero resistance. The equations of motion relate the time dependence 
of the Goldstone boson field to the potential, 



(x) = -Vix). 



(53) 



The electric current is related to the gradient of the Goldstone boson field. Equ. (|53|l shows that the time dependence 
of the current is proportional to the gradient of the potential. In order to have a static current the gradient of the 
potential has to be constant throughout the sample, and the resistance is zero. 

In order to study the properties of a superconductor in more detail we have to specify Ls- For this purpose we 
assume that the system is time-independent, that the spatial gradients are small, and that the order parameter is 
small. In this case we can write 



Ls= d-" 



(54) 



where uih and g are unknown parameters that depend on the temperature. Equ. (|54|l is known as the Landau- 
Ginzburg effective action. Strictly speaking, the assumption that the order parameter is small can only be justified 
in the vicinity of a second order phase transition. Nevertheless, the Landau-Ginzburg description is instructive even 
in the regime where t = {T — Tc)/Tc is not small. It is useful to decompose $ — pexp{2ie(p). For constant fields the 
effective potential. 



Vip) 



1 2 2,1 4 



(55) 



is independent of (f>. The minimum is at Pq — nijj / g and the energy density at the minimum is given by £ = ""itt-h I (^ff) • 
This shows that the two parameters m/f and g can be related to the expectation value of $ and the condensation 
energy. We also observe that the phase transition is characterized by m^f (T^) = 0. 
In terms of (/) and p the Landau-Ginzburg action is given by 
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(56) 



The equations of motion for A and p are given by 

V X B 



4eV^ (V0 - a) 



2 2 



99 



4eV A . 



(57) 
(58) 



Equ. (EU imphes that V^B = -'^e^p^B. This means that an external magnetic field B decays over a characteristic 
distance A = l/(2ep). Equ. (|58|) gives V^p — —'rnjjp+ .... As a consequence, variations in the order parameter relax 
over a length scale given by ^ = l/mn- The two parameters A and ^ are known as the penetration depth and the 
coherence length. 

The relative size of A and ^ has important consequences for the properties of superconductors. In a type II 
superconductor ^ < A. In this case magnetic flux can penetrate the system in the form of vortex lines. At the core 
of a vortex the order parameter vanishes, p = 0. In a type II material the core is much smaller than the region over 
which the magnetic field goes to zero. The magnetic flux is given by 



A-dl ^ f V(t)-dl^ 
dA JdA e 



(59) 
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and quantized in units of -kTi/ e. In a type II superconductor magnetic vortices repel each other and form a regular 
lattice known as the Abrikosov lattice. In a type I material, on the other hand, vortices are not stable and magnetic 
fields can only penetrate the sample if superconductivity is destroyed. 

The Landau-Ginzburg description shows that there is no qualitative difference between superconductivity and Bose 
condensation of charged bosons. Indeed, we may think of superconductivity as Bose condensation of Cooper pairs. 
While this is qualitatively correct, there is an important quantitative difference between a BCS superconductor and 
a dilute Bose condensate of composite bosons. In a BCS superconductor the coherence length ^, which is a measure 
of the size of the Cooper pairs, is much larger than the average inter-particle spacing p^^. Also, the pair correlation 
essentially disappears above the critical temperature. In a dilute Bose condensate, on the other hand, the size of the 
bosons is much smaller than the typical distance between them. The bosons are tightly bound and do not dissolve at 
Tc- Nevertheless, since there is no qualitative difference between Bose condensation and BCS superconductivity we 
expect to find systems that show a crossover from one kind of behavior to the other. We will discuss an example in 
Sect. lVIIIDI 



V. QCD AND SYMMETRIES 



Before we discuss QCD at finite baryon density we would like to provide a quick reminder on QCD and the 
symmetries of QCD. The elementary degrees of freedom are quark fields ■0° ^ and gluons A'^^. Here, a is color index 
that transforms in the fundamental representation for fermions and in the adjoint representation for gluons. Also, / 
labels the quark flavors u, d, s, c, 6, t. In practice, we will focus on the three light flavors up, down and strange. The 
QCD lagrangian is 

Nf _ 

C = J2 M^P - ^f)^f - -^Gl^Gl,, (60) 
/ 

where the field strength tensor is defined by 

G^,. = d^Al - d,A1 + gr'-^AlAl, (61) 
and the covariant derivative acting on quark fields is 



iipi,^Y\^d^+9Al — j^. (62) 

QCD has a number of remarkable properties. Most remarkably, even though QCD accounts for the rich phenomenology 
of hadronic and nuclear physics, it is an essentially parameter free theory. To first approximation, the masses of the 
light quarks m, d, s are too small to be important, while the masses of the heavy quarks c, 5, t are too heavy. If we set 
the masses of the light quarks to zero and take the masses of the heavy quarks to be infinite then the only parameter 
in the QCD lagrangian is the coupling constant, g. Once quantum corrections are taken into account g becomes a 
function of the scale at which it is measured. If the scale is large then the coupling is small, but in the infrared the 
coupling becomes large. This is the famous phenomenon of asymptotic freedom. Since the coupling depends on the 
scale the dimensionless parameter g is traded for a dimensionful scale parameter Aqcd- In essence, Aqcd is the scale 
at which the coupling becomes large. 

Since Aqcd is the only dimensionful quantity in QCD (m^ = 0) it is not really a parameter of QCD, but reflects 
our choice of units. In standard units, Aqcd — 200 MeV ~ 1 fm"""^. Note that hadrons indeed have sizes r ~ Aq]-,^. 
However, we should also note that in practice the perturbative expansion in g breaks down at scales r ^ ^xSB ^ 
0.2fm«AQj,^. 

Another important feature of the QCD lagrangian are its symmetries. First of all, the lagrangian is invariant under 
local gauge transformations U{x) £ SU{3)c 

ip{x) ^ U{x)ipix), A^{x) U{x)A^U''{x) + iU{x)d^,U\x), (63) 

where A^ = ^J^(A°/2). In the QCD ground state at zero temperature and density the local color symmetry is confined. 
This implies that all excitations are singlets under the gauge group. 

The dynamics of QCD is completely independent of fiavor. This implies that if the masses of the quarks are equal, 
ruu — rrid — rris , then the theory is invariant under arbitrary flavor rotations of the quark fields 

V'/ ^ Vfg^g, (64) 
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where V € SU{3). This is the weh known flavor (isospin) symmetry of the strong interactions. If the quark masses 
are not just equal, but equal to zero, then the flavor symmetry is enlarged. This can be seen by defining left and 
right-handed fields 

ijL,R = l{l±l5)i^. (65) 

In terms of L/R fields the fermionic lagrangian is 

£ = i)L{iiP)i)L + i)R{iip)i}R + i)LMi,R + i^RM^jL, (66) 

where M = diag(m„, m^, to^). We observe that if quarks are massless, m„ = = = 0, then there is no 
coupling between left and right handed fields. As a consequence, the lagrangian is invariant under independent flavor 
transformations of the left and right handed fields. 

tpLJ Lfglj;L,g, IpRJ Rfgi^R,g, (67) 

where {L,R) e SU{3)l x SU{3)r. In the real world, of course, the masses of the up, down and strange quarks are 
not zero. Nevertheless, since m„,md <^ms < Aqcd QCD has an approximate chiral symmetry. 

In the QCD ground state at zero temperature and density the flavor symmetry is realized, but the chiral symmetry 
is spontaneously broken by a quark-anti-quark condensate {^jl^jr + iJjr.iJjl)- As a result, the observed hadrons can be 
approximately assigned to representations of the SU{3)v flavor group, but not to representations of SU (3) l x SU{3)ii. 
Nevertheless, chiral symmetry has important implications for the dynamics of QCD at low energy. Goldstone's 
theorem implies that the breaking of SU{3)l x SU{3)u SU{3)v is associated with the appearance of an octet 
of (approximately) massless pseudoscalar Goldstone bosons. Chiral symmetry places important restrictions on the 
interaction of the Goldstone bosons. These constraints are obtained most easily from the low energy effective chiral 
lagrangian. At leading order we have 



^Tr [df.^d"!:''] + STr(MEt) + h.c. +..., (68) 



where S = exp{i(f>°- X'^ / fj^) is the chiral field, /,r is the pion decay constant and M is the mass matrix. Expanding E 
in powers of the pion, kaon and eta fields we can derive the leading order chiral perturbation theory results for 
Goldstone boson scattering and the coupling of Goldstone bosons to external fields. Higher order corrections originate 
from loops and higher order terms in the effective lagrangian. 

Finally, we observe that the QCD lagrangian has two U{1) symmetries, 

C/(1)b: V'L^e'^^L, ^Pr^c^^^Pr (69) 

C/(1)a: ^L^e'«Vi, ijR^e-^"^R. (70) 

The U{1)b symmetry is exact even if the quarks are not massless. Superficially, it appears that the U{1)a symmetry 
is explicitly broken by the quark masses and spontaneously broken by the quark condensate. However, there is no 
Goldstone boson associated with spontaneous U{1)a breaking. The reason is that at the quantum level the U{1)a 
symmetry is broken by an anomaly. The divergence of the U{1)a current is given by 

where G^j^ — is the dual field strength tensor. 



VI. QCD AT FINITE DENSITY 

In the real world the quark masses are not equal and the only exact global symmetries of QCD are the U{l)f 
flavor symmetries associated with the conservation of the number of up, down, and strange quarks. If we take into 
account the weak interactions then flavor is no longer conserved and the only exact symmetries are the U{1)b of 
baryon number and the U{1)q of electric charge. 

In the following we study hadronic matter at non-zero baryon density. We will mostly focus on systems at non-zero 
baryon chemical potential but zero electron U{1)q chemical potential. We should note that in the context of neutron 
stars we are interested in situations when the electric charge, but not necessarily the electron chemical potential, is 
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FIG. 6: Naive phase diagram of hadronic matter as a function of the baryon chemical potential and temperature. 



zero. We will comment on the consequences of electric charge neutrality below. Also, if the system is in equilibrium 
with respect to strong, but not to weak interactions, then non-zero flavor chemical potentials may come into play. 
The partition function of QCD at non-zero baryon chemical potential is given by 

where i labels all quantum states of the system, Ei and are the energy and baryon number of the state i. If the 
temperature and chemical potential are both zero then only the ground state contributes to the partition function. 
All other states give contributions that are exponentially small if the volume of the system is taken to infinity. In 
QCD there is a massgap for states that carry baryon number. As a a consequence there is an onset chemical potential 

llonset^^VCi{EilN.i), (73) 

% 

such that the partition function is independent of /i for /i < [lonset- For \x > [lonset the baryon density is non-zero. 
If the chemical potential is just above the onset chemical potential we can describe QCD, to first approximation, as 
a dilute gas of non- interacting nucleons. In this approximation ii onset = m^. Of course, the interaction between 
nucleons is essential. Without it, we would not have stable nuclei. As a consequence, nuclear matter is self-bound 
and the energy per baryon in the ground state is given by 

=i -15MeV. (74) 

The onset transition is a first order transition at which the baryon density jumps from zero to nuclear matter saturation 
density, po ~ 0.14fm~"^. The first order transition continues into the finite temperature plane and ends at a critical 
endpoint at T = ~ 10 MeV, see Fig. 

Nuclear matter is a complicated many-body system and, unlike the situation at zero density and finite temperature, 
there is also no information from numerical simulations on the lattice. This is related to the so-called 'sign problem'. 
At non-zero chemical potential the euclidean fermion determinant is complex and standard Monte-Carlo techniques 
based on importance sampling fail. Recently, some progress has been made in simulating QCD for small /i and T ~ Tc 
[2^ I271 1281] ■ but the regime of small temperature remains inaccessible. As a consequence of the sign problem, there 
are also essentially no general results concerning the structure of the ground state. While the theorems of Vafa and 
Witten [2^ Isof rule out spontaneous breaking of parity or flavor at T 7^ and 11 = 0, there are no theorems of this 
type at non-zero baryon density. 

However, if the density is very much larger than nuclear matter saturation density, po,we expect the problem to 
simplify. In this regime it is natural to use a system of non-interacting quarks as a starting point | 3l| . The low energy 
degrees of freedom are quark excitations and holes in the vicinity of the Fermi surface. Since the Fermi momentum is 
large, asymptotic freedom implies that the interaction between quasi-particles is weak. As a consequence, the naive 
expectation is that chiral symmetry is restored and quarks and gluons are deconfined. It seems natural to assume that 
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the quark liquid at high baryon density is continuously connected to the quark-gluon plasma at high temperature. 
These naive expectations are summarized in the phase diagram shown in Fig. 

Corrections to the non-interacting quark liquid can be studied in perturbation theory. The thermodynamic potential 
is given by [s^, 



47r2 



where G = Gq — 0.536iV/ + Nf \nNf, Gq = 10.374 ± 0.13. Here, /i is the chemical potential for quark number. This 
convention is more natural in the context of perturbative QCD and we will use it for the remainder of these lectures. 
Note that perturbative corrections reduce the pressure of the quark phase. At least qualitatively, this is agreement 
with the idea that at very low density the pressure of the hadron phase is bigger than the pressure of the quark phase. 



VII. COLOR SUPERCONDUCTIVITY 



There are two problems with the perturbative expansion equ. (|75|) . One problem is related to the fact that while the 
electric gluon interaction is screened by the mechanism discussed in Sect. HlBl there is no screening of magnetic gluon 
exchanges. This not only implies that the mag netic sector of the theory becomes non-perturbative, it also causes the 
Fermi liquid description to break down 1^J35| . The correction to the fermion self energy near the Fermi surface due 



to magnetic gluon exchanges is |36l I37ll38l l3 



This correction invalidates the Fermi liquid description for energies po ~ fj,exp{—l/g^). But even before this phe- 
nomenon becomes important there is another effect that will invalidate the Fermi liquid picture. In Sect. we 
showed that the BCS instability will lead to pair condensation whenever there is an attractive fermion-fermion interac- 
tion. At very large density, the attraction is provided by one-gluon exchange between quarks in a color anti-symmetric 
3 state. High density quark matter is therefore expected to behave as a color superconductor |4nL lill 142. l43l | . 
Color superconductivity is described by a pair condensate of the form 

$ = (^^CFdActfV')- (77) 

Here, C is the charge conjugation matrix, and Tu, Xc,tf are Dirac, color, and flavor matrices. Except in the case 
of only two colors, the order parameter cannot be a color singlet. Color superconductivity is therefore characterized 
by the breakdown of color gauge invariance. This statement has to be interpreted in the sense of Sect. IIVUI Gluons 
acquire a mass due to the (Meissner- Anderson) Higgs mechanism. 

A rough estimate of the critical density for the transition from chiral symmetry breaking to color superconductivity, 
the superconducting gap and the transition temperature is provided by schematic four-fermion models ,44..45j . Typical 
models are based on the instanton interaction 

C^Gi{{i,T-yjf + {^^,r^i,r}, (78) 

or a schematic one-gluon exchange interaction 

X 2 



C = GoGE [i'l^Y'^ j ■ (79) 

Here — (r, «) is an isospin matrix and A" are the color Gell-Mann matrices. The strength of the four-fermion 
interaction is typically tuned to reproduce the magnitude of the chiral condensate and the pion decay constant at zero 
temperature and density. In the mean field approximation the effective quark mass associated with chiral symmetry 
breaking is determined by a gap equation of the type 

= Gm I It \, , ^'^^ (1 - nF{E,)) , (80) 



where Gm is the effective coupling in the quark-anti-quark channel and A is a cutoff. Both the instanton interaction 
and the one-gluon exchange interaction are attractive in the color anti-triplet scalar diquark channel e°'^''{4>''Cj5^'^). 
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FIG. 7: First revision of the phase diagram of hadronic matter. This figure shows the phase diagram of strongly interacting 
matter obtained from a mean field treatment of chiral symmetry breaking and color superconductivity in QCD with two flavors, 
see e.g. [Hof. 

A pure one-gluon exchange interaction leads to a degeneracy between scalar and pseudoscalar diquark condensation, 
but instantons are repulsive in the pseudoscalar diquark channel. The gap equation in the scalar diquark channel is 



where we have neglected terms that do not have a singularity on the Fermi surface, \p\ ^ pp- In the case of a four- 
fermion interaction with the quantum numbers of one-gluon exchange Go = Gm/{Nc — 1). The same result holds 
for instanton effects. In order to determine the correct ground state we have to compare the condensation energy in 
the chiral symmetry broken and diquark condensed phases. We have £ ~ f^-^'^Q the {qq) condensed phase and 
£ ^ p|,A^/(27r^) in the {qq) condensed phase. 

At zero temperature and density both equs. H8l)|l and H81|l only have non-trivial solutions if the coupling exceeds 
a critical value. Since Cm > Gd have Mq > A and the energetically preferred solution corresponds to chiral 
symmetry breaking. If the density increases Pauli-Blocking in equ. H80|l becomes important and the effective quark 
mass decreases. The diquark gap equation behaves very differently. Equ. H81|l has an infrared singularity on the Fermi 
surface, p = pp, and this singularity is multiplied by a finite density of states, N — p|,/(27r)^. As a consequence, there 
is a non-trivial solution even if the coupling is weak. The gap grows with density until the Fermi momentum becomes 
on the order of the cutoff. For realistic values of the parameters we find a first order transition for remarkably small 
values of the quark chemical potential, fiQ ~ 300 MeV. The gap in the diquark condensed phase is A ~ 100 MeV and 
the critical temperature is Tc ~ 50 MeV. 

In the same model the finite temperature phase transition at zero baryon density is found to be of second order. 
This result is in agreement with universality arguments 46] and lattice results. If the transition at finite density 
and zero temperature is indeed of first order then the first order transition at zero baryon density has to end in a 
tri-critical point 113, 113) 13 • The tri-critical point is quite remarkable, because it remains a true critical point, 
even if the quark masses are not zero. A non-zero quark mass turns the second order T ^ transition into a smooth 
crossover, but the first order /i ^ transition persists. While it is hard to predict where exactly the tri-critical point 
is located in the phase diagram it may well be possible to settle the question experimentally. Heavy ion collisions at 
relativistic energies produce matter under the right conditions and experimental signatures of the tri-critical point 
have been suggested in i52] . 

A schematic phase diagram is shown in Fig. [7| We should emphasize that this phase diagram is based on simplified 
models and that there is no proof that the transition from nuclear matter to quark matter along the T = line occurs 
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FIG. 8: Dyson-Schwinger (gap) equation in QCD at finite density. The square denotes an anomalous self energy (gap) insertion, 
and the tria ngle denotes a gluon self energy insertion. At leading order, quark self energy insertions or vertex corrections are 
not required [23 • 

via a single first order transition. Chiral symmetry breaking and color superconductivity represent two competing 
forms of order, and it seems unlikely that the two phases are separated by a second order transition. However, since 
color superconductivity modifies the spectrum near the Fermi surface, whereas chiral symmetry breaking operates 
near the surface of the Dirac sea, it is not clear that the two phases cannot coexist. Indeed, there are models in which 
a phase coexistence region appears [s^ . 



VIII. PHASE STRUCTURE IN WEAK COUPLING 



A. QCD with two flavors 



In this section we shall discuss how to use weak coupling methods in order to explore the phases of dense quark 
matter. We begin with what is usually considered to be the simplest case, quark matter with two degenerate flavors, up 
and down. Renormalization group arguments suggest 54, 55], and explicit calculations show ,39, 56], that whenever 
possible quark pairs condense in an s-wave. This means that the spin wave function of the pair is anti-symmetric. 
Since the color wave function is also anti-symmetric, the Pauli principle requires the flavor wave function to be 
anti-symmetric too. This essentially determines the structure of the order parameter 0, 

This order parameter breaks the color SU{3) —^ SU(2) and leads to a gap for up and down quarks with two out of 
the three colors. Chiral and isospin symmetry remain unbroken. 

We can calculate the magnitude of the gap and the condensation energy using weak coupling methods. In weak 
coupling the gap is determined by ladder diagrams with the one gluon exchange interaction. These diagrams can be 
summed using the gap equation [53, |5^ |5^ [Sfl 0] 

9^ r . r . . ( |-icos0 



A(.4) -^Jdl^J ^cos. [j^i^^^f^^ (83) 



l-cose + F/{2^?)) + A(g4)2 ' 

Here, A(p4) is the frequency dependent gap, g is the QCD coupling constant and G and F are the self energies 
of magnetic and electric gluons. This gap equation is very similar to the BCS gap equation equ. (|81|l obtained in 
four-fermion models. The terms in the curly brackets arise from the magnetic and electric components of the gluon 
propagator. The numerators are the on-shell matrix elements Mu^qq = [uh{pi)^i,QUh{p?,)][uh{p2)li,aUh{p4)] for the 
scattering of back-to-back fermions on the Fermi surface. The scattering angle is cos = pi ■ p^. In the case of a spin 
zero order parameter, the helicity h of all fermions is the same, see for more detail. 

The main difference between equ. (|83|) and the BCS gap equation (|81|) is that because the gluon is massless, the 
gap equation contains a coUinear cos 0^1 divergence. In a dense medium the coUinear divergence is regularized by 
the gluon self energy. For (f — > and to leading order in perturbation theory we have 

F = 2m2, G^-m^^, (84) 

2 Jg]' ^ ' 



with = N f fj? / {^tt"^) . In the electric part, m|, = 2m? is the familiar Debye screening mass. In the magnetic 
part, there is no screening of static modes, but non-static modes are modes are dynamically screened due to Landau 
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damping. Equ. H84|l is, up to an overall degeneracy factor, exactly equal to the result obtained in Sect. Ill Bl The 
only difference is that in a relativistic theory the role of the tadpole graph in Fig. |^ is played by the contribution 
of negative energy states in the particle-hole graph Fig. Et- We refer the reader to p2i lo l64l | for a more complete 
discussion of quasi-particle properties in a dense quark liquid. 

For small energies dynamic screening of magnetic modes is much weaker than Debye screening of electric modes. 
As a consequence, perturbative color superconductivity is dominated by magnetic gluon exchanges. Using equ. H84|l 
we can perform the angular integral in equ. (|83|) and find 

A(P4) -4^fd,. log f -j^] -r^^, (85) 

with b = 2567r^(2/Af/)^/^g^^. We can now see why it was important to keep the frequency dependence of the gap. 
Because the collinear divergence is regulated by dynamic screening, the gap equation depends on p4 even if the 
frequency is small. We can also see that the gap scales as exp(— c/g). The collinear divergence leads to a gap equation 
with a double-log behavior. Qualitatively 



1 
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log (86) 



from which we conclude that A exp(—c/g). The approximation equ. (|86|l is not sufficiently accurate to determine 
the correct value of the constant c. A more detailed analysis shows that the gap on the Fermi surface is given by 



(87) 



The factor 6q is related to non-Fermi liquid effects, see equ. (|76ll . Note that since A ^ exp(— 1/q) non-Fermi liquid 
effects are indeed sub-leading. In perturbation theory 6q = cxp(— (tt^ -|- 4)(iVc — 1)/16) |6j |63- The condensation 
energy is given by 



.2 / 



where Nd = 4 is the number of condensed species. The critical temperature is Tc/Aq — j-K ~ 0.56, as in standard 
BCS theory. For chemical potentials /i < 1 GeV, the coupling constant is not small and the applicability of perturba- 
tion theory is in doubt. If we ignore this problem and extrapolate the perturbative calculation to densities p ~ 5po we 
find gaps A ~ 100 MeV. This result is in surprisingly good agreement with the estimates from Nambu-Jona-Lasinio 
models discussed in Sect. IVIII 

We note that the 2SC phase defined by equ. has two gapless fermions and an unbroken SU{2) gauge group. 
The gapless fermions are singlets under the unbroken SU(2). As a consequence, we expect the SU(2) gauge group 
to become non-perturbative. An estimate of the SU{2) confinement scale was given in [6g. We also note that even 
though the Copper pairs carry electric charge the U{1) of electromagnetism is not broken. The generator of this 
symmetry is a linear combination of the original electric charge operator and the diagonal color charges. Under this 
symmetry the gapless fermions carry the charges of the proton and neutron. Possible pairing between the gapless 
fermions was discussed in 1441671. 



B. QCD with three flavors: Color-Flavor-Locking 



If quark matter is formed at densities several times nuclear matter density we expect the quark chemical potential to 
be larger than the strange quark mass. We therefore have to determine the structure of the superfluid order parameter 
for three quark flavors. We begin with the idealized situation of three degenerate flavors. From the arguments given 
in the last section we expect the order parameter to be color and flavor anti-symmetric matrix of the form 

<i>^'' = (CC75V'-)- (89) 
In order to determine the precise structure of this matrix we have to extremize grand canonical potential. We find 

mil 



A^ = AAiSIS] - S'^S^) + A5(<5^5^ + S'^S]), 



(90) 
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FIG. 9: Conjectured phase diagram of Nf = 3 hadronic matter in the limit of exact flavor symmetry. 



which describes the color-flavor locked (CFL) phase proposed in |70|- In the weak coupling limit A5 ^ A^i and 
= 2~^/^Ao where Ag is the gap in the 2SC phase, equ. IjHTII In the CFL phase both color and flavor symmetry 
are completely broken. There are eight combinations of color and flavor symmetries that generate unbroken global 
symmetries. The unbroken symmetries are 

rL,^ - (U^rU^.^h, rn,. - (.U*rU.,i:'k^, (91) 

for U e SU{3)v- The symmetry breaking pattern is 

SU{3)l-x SU{3)rxU{1)v ^ SU{3)v (92) 

We observe that color-flavor-locking implies that chiral symmetry is broken. The mechanism for chiral symmetry 
breaking is quite unusual. The primary order parameter {ij^l-C IS.'^^^^-) = —{ip'^-CAf-'ip'^j^-) involves no coupHng 
between left and right handed fermions. In the CFL phase both left and right handed flavor are locked to color, and 
because of the vectorial coupling of the gluon left handed flavor is effectively locked to right handed flavor. Chiral 
symmetry breaking also implies that (tpip) has a non-zero expectation value. We shall compute the quark condensate 
in Sect. lIXEl In the CFL phase {•ipip)'^ ^ {{ipip)'^). Another measure of chiral symmetry breaking is provided by the 
pion decay constant. In Sect. HXDl we will show that in the weak coupling limit is proportional to the density of 
states on the Fermi surface. 

The symmetry breaking pattern SU{3)l x SU{3)u SU{3)v in the CFL phase is identical to the symmetry 
breaking pattern in QCD at low density. The spectrum of excitations in the color-flavor-locked (CFL) phase also 
looks remarkably like the spectrum of QCD at low density jT^]. The excitations can be classified according to their 
quantum numbers under the unbroken SU(3), and by their electric charge. The modified charge operator that 
generates a true symmetry of the CFL phase is given by a linear combination of the original charge operator Qem and 
the color hypercharge operator Q = diag(— 2/3, —2/3, 1/3). Also, baryon number is only broken modulo 2/3, which 
means that one can still distinguish baryons from mesons. We find that the CFL phase contains an octet of Goldstone 
bosons associated with chiral symmetry breaking, an octet of vector mesons, an octet and a singlet of baryons, and a 
singlet Goldstone boson related to superfluidity. All of these states have integer charges. 

With the exception of the U{1) Goldstone boson, these states exactly match the quantum numbers of the lowest 
lying multiplets in QCD at low density. In addition to that, the presence of the U{1) Goldstone boson can also be 
understood. The U{1) order parameter is {{uds) (uds)) . This order parameter has the quantum numbers of a 0"*" AA 
pair condensate. In Nf — 3 QCD, this is the most symmetric two nucleon channel, and a very likely candidate for 
superfluidity in nuclear matter at low to moderate density. We conclude that in QCD with three degenerate light 
flavors, there is no fundamental difference between the high and low density phases. This implies that a low density 
hyper-nuclear phase and the high density quark phase might be continuously connected, without an intervening phase 
transition. A conjectured phase diagram is shown in Fig. |51 
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FIG. 10: Schematic plot of the properties of hadronic modes in Nf = 3 QCD in the hmit of exact flavor symmetries. Hadronic 
parameters can be determined using chiral perturbation theory at low density, and using perturbative QCD at high density. 
The high density behavior is studied in Sect. IIX'HI 

An important consistency check for the structure of the phase diagram is provided by anomaly matching arguments. 
Anomaly matching expresses the requirement that the flavor anomalies of the microscopic theory can be represented 
in the effective theory for the low energy degrees of freedom [f^, IZ^ ■ It was shown in 0, that this requirement 
also applies to gauge theories at finite baryon density. We observe that color-flavor-locking realizes the standard 
Goldstone boson option for anomaly matching in Nf = 3 QCD, whereas the 2SC phase corresponds to the massless 
proton and neutron option in Nf — 2 QCD. 

We also note that the CFL phase provides a weak coupling realization of a phase with chiral symmetry breaking 
and a mass gap. This means that the CFL phase offers the opportunity to study many of the 'hard' problems of 
non-perturbative QCD in a perturbative setting. For example, we can compute hadronic parameters at high baryon 
density and try to extrapolate the results to low density, see Fig. ^| 

C. A^//2,3 

Color-flavor locking can be generalized to QCD with more than three flavors IgS] . For all Nf > 3 the high density 
phase is fully gapped. Also, at least part of the chiral SU{Nf) x SU{Nf) symmetry is broken for all Nf > 3, but only 
in the case Nf = 3 do we find the T = /i = pattern of chiral symmetry breaking, SU{Nf)]^ x SU{Nf)fi SU{Nf)Y. 

While the case Nf > 3 is mostly of academic interest, the phase structure of Nf = 1 QCD is possibly relevant to real 
QCD. Because of mass effects or non-zero electron chemical potentials the Fermi surface of strange and non-strange 
quarks may get pushed too far apart for strange-non-strange pairing to occur, see Sect. lIX As a consequence, there 
may be regions in the phase diagram where (ss) pairing occurs. In the case of single flavor pairing the order parameter 
is flavor-symmetric and the Cooper pairs carry non-zero angular momentum. The simplest order parameters are of 
the form 

Si = (e"''"?A''C7V'"), $^ = (e''''"V''C#"). (93) 

The corresponding gaps can be determined using the methods introduced in section IVlII Al We find A($i 2) = 
exp(— 3ci^2) Aq where ci = —1.5, C2 — —2 and Aq is the spin zero gap [s^jH^. While the natural scale of the s-wave 
gap is Aq — 100 MeV, the p-wave gap is expected to be less than 1 MeV. 

The spin one order parameter equ. H93|) is a color-spin matrix. This opens the possibility that color and spin 
degrees become entangled, similar to the color-flavor-locked phase or the B-phase of liquid '^He. The corresponding 
order parameter is 

^CSL = Sfie^'^tP'O {cos{(3)q^ + MPH) r), (94) 

where the angle f3 determines the mixing between the two types of condensates shown in equ. (|93|l . A weak coupling 
analysis of the effective potential shows that the color-spin-locked phase equ. H94|l is favored over the "polar" phase 
equ. I|93|l [s^ . The value of /? depends sensitively on the interaction and the mass of the quark. In the color-spin- 
locked phase color and rotational invariance are broken, but a diagonal S0{3) survives. As a consequence, the gap is 
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isotropic. Color-spin-locking also leads to an unusual spectrum of quasi-particles. In the Fermi liquid phase there is 
a left and right handed color triplet of quarks. In the CSL phase we find a spin 3/2 quartet and a spin 1/2 doublet 
of the unbroken 50(3) symmetry. The CSL phase 'knows' that one-flavor nuclear matter consists of spin 3/2 delta 
baryons. 



D. N^ = 2 



QCD with Nc~ 2 colors is an interesting model system. The interest in this theory derives from the fact that the 
determinant of the euclidean Dirac operator in QCD with Nc = 2 remains real even if the baryon chemical potential 
is non-zero. This means that two-color QCD with an even number of flavors can be studied on the lattice using 
standard techniques. In addition to that, there is theoretical control not only in the regime of large density, but also 
in the regime of small density. 

For simplicity we will concentrate on the case of Nf = 2 flavors. SU{2) gauge theory has a meson spectrum which 
is very similar to three-color QCD. Baryons, on the other hand, are bosons rather than fermions and their spectrum 
is very different as compared to iVc = 3 QCD. Two-color QCD is also characterized by an enlarged chiral symmetry. 
We can write 



i'L 



(95) 



where (T2, T2 are anti-symmetric color and flavor SU(2) matrices. Two-color QCD is not only invariant under SU (2)^ x 
SU(2)ji transformations acting on the upper and lower components of separately, but under the full S'C/(4) group 
B^i2fi.X7]- The SU{4:) chiral symmetry mixes the quark-anti-quark condensate {'ip'^ip'^) with the diquark condensate 

At zero temperature and density, and in the presence of a small quark mass, the chiral SU{^) symmetry is broken to 
Sp{A) by a quark-anti-quark condensate {{^"■■ip"') . There are 5 Goldstone bosons, three pions tt , the scalar diquark S and 
the scalar anti-diquark S. If we turn on a baryon chemical potential then the scalar diquark will Bose condense if the 
chemical potential exceeds the mass of the diquark. Since the scalar diquark is a Goldstone boson, this phenomenon 
can be studied using the chiral effective lagrangian of A^c = 2 QCD. The effective lagrangian is given by [t^F^ 



C = ^Tr [D^SD^S^] + TOCTr(AfEt) + /^.^ 



(96) 



The chiral field E parametrizes the coset SU{4)/ Sp{A). E is an anti-symmetric unitary matrix. 
S — > UT,U^ under the chiral SU (4) symmetry. The covariant derivative is defined as 



It transforms as 



(97) 



The matrices M and B are determined by the transformation properties of the mass term mipijj and the chemical 
potential term jiip'^'il! under the chiral SU{4) symmetry. We have 



M 



1 

-1 



B 



1 

-1 



(98) 



We can determine the ground state by minimizing the potential 

J TT 



F(E) 



Tr [ES^Etfi] + TOCTr(AfE^) + h.c. 



(99) 



For zero density and finite quark mass the minimum is E^^ = M. In this state S'C/(4) is broken to Sp{A). Goldstone 
bosons are described by fluctuations E = UY^qqlJ^ with U = cxp{i(j)''X'' / e S'C/(4)/S'p(4). Here, are the S'[/(4) 
generators that act non-trivially on E^g. As mentioned above, there a 5 Goldstone modes. 
The first term in equ. (|99|l is minimized by 



Sgq — 



<J2 

a2 



(100) 
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For non-zero chemical we expect the minimum to be of the form Eg = cos(a) -I- T,qqS'm{a). Substituting this 
ansatz into equ. H99|l we find a = for /i < 171^^/2 and cos(a) = m^/(4/i)^. Differentiating the effective potential with 
respect to the chemical potential we find the baryon density 



PB = 8/2^1- . (101) 




For small /i — m^/2+(5/i this result is exactly of the same form as equ. H27(l . This implies that the physical phenomenon 
is indeed Bose condensation of scalar diquarks interacting through a short range repulsive interaction. 

If the density becomes large, /i m,p ^ m^, the effective lagrangian description breaks down. On the other hand, if 
/i 3> ^QCD we expect to find a perturbative BCS superfluid of diquark pairs. The gap and the superfluid condensate 
can be computed using the methods discussed in Sect. IVIII Al We find 

A = bl2-K%ijLg-'^ exp ' (102) 

with 6q = exp(— (tt^ -|- 4)/16). Numerical studies on the lattice can be used in order to verify the limiting behavior at 
small and large density, and to study the Bose condensation/BCS crossover, see Q| and references therein. Similar 
studies have also been performed in the instanton liquid model |8C| . There are many other interesting questions 
that can be studied in two-color QCD. It was suggested, for example, that vector mesons (diquarks) may condense 
if the chemical potential is on the order of the vector meson mass 81, 82J. Other interesting questions concern the 
structure of the phase diagram at non-zero temperature [s^, and the nature of the deconfinement transition. It was 
also pointed out that the behavior of the rj' mass in two-color QCD can be used in order to study the mechanism 
of U{1)a breaking in QCD 84]. Finally, we should mention that there are some other gauge theories in which the 
euclidean fermion determinant is positive even if the chemical potential is non-zero. These theories include QCD at 
finite isospin chemical potential and QCD with a non-zero density of quarks in the adjoint representation of color 
[tsL Issl ISq . It is amusing that all of these theories have Goldstone bosons that carry the conserved charge. As a 
consequence, the low density state of these theories is a dilute Goldstone boson condensate, similar to the Nc ^ 2 
diquark condensate studied in this section. 

E. Nc^ oo 

In the large Nc limit quark-quark scattering is suppressed as compared to quark-quark-hole scattering. The color 
factors in the two channels are 

+ 1 A^ - 1 

c=^ (..), (,.), (103) 

suggesting that particle-particle pairing, and superconductivity, is disfavored at large Nc- As a consequence, other 
forms of pairing may take place. Particle-hole scattering is not suppressed, but particle-hole pairing can only take 
place over a small part of the Fermi surface. The order parameter for particle-hole pairing is |87j 

{i>{x)ij{y)) = exp(zp • (f + y))Y.{x - y), (104) 

where \p\ — pp is a vector on the Fermi surface. This state describes a chiral density wave. It resembles a charge 
or spin density wave in quasi-one-dimensional condensed matter systems js^ . QCD, of course, is not quasi-one- 
dimensional but at large Nc screening due to fermions is weak and the perturbative one-gluon exchange interaction 
is very strongly dominated by coUinear scattering. The problem was studied in the weak coupling approximation in 
[89j |. It was found that the transition from color superconductivity to chiral density waves requires very large values 
of Nc > 1000. On the other hand, if the coupling is strong and the density is not too large, then the chiral density 
wave state may compete with color superconductivity even for three colors [90|. 

IX. THE ROLE OF THE STRANGE QUARK MASS 
A. BCS theory: toy model 



At baryon densities relevant to astrophysical objects d istor t ions of the pure CFL state due to non-zero quark masses 
cannot be neglected [llilililiEilillgiil [loS EoJ- The most important eflFect of a non-zero strange 
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FIG. 11: Light quark gap A^d and strange-non-strange gap Aua as a function of the strange quark mass for p = 0.5 GeV. 
This figure shows the result of a mean field analysis based on a schematic BCS interaction ^J. The interaction strength was 
(arbitrarily) adjusted to A„tj(ms=0) — 25 MeV. 



quark mass is that the light and strange Fermi surfaces will no longer be of equal size. When the mismatch is much 
smaller than the gap one calculates assuming degenerate quarks, we might expect that it has very little consequence, 
since at this level the original particle and hole states near the Fermi surface are mixed up anyway. On the other hand, 
when the mismatch is much larger than the nominal gap, we might expect that the ordering one would obtain for 
degenerate quarks is disrupted, and that to a first approximation one can treat the light and heavy quark dynamics 
separately. 

We can see this in a more quantitative fashion by studying a schematic gap equation that describes the spin singlet 
pairing of two fermions with different masses. In a basis of particles of the first kind and holes of the second the 
quadratic part of the action is 



S 



(2^ 



(V^(l) ^(2) 



A 



A* 



Po 



^(1) 



(105) 



Here, Cp ^ = E'p'^ — and E'^''^ = (p^ -I- 2)^^^ where mi_2 are the masses of particle one and two. The particle and 
hole propagators are determined by the inverse of the matrix equ. H1U5|) . The off-diagonal (anomalous) propagator is 

G2I = 7 ^ . o^ —■ (106) 



{p,-el){po + el)-^^■ 

We study the effect of a zero range interaction G'(V'iO'2V'2)('0i<''2V'2)- The pairing is described by the gap equation 

d-^p A 



A = G 



(2^)4 (po + i? + ifcgn(po))2 - e2 - A2 ■ 



(107) 



Here, we have introduced ep — Ep — fi — {ej, + e^)/2 and R — {tp — ep)/2. In practice, we are interested in pairing 
between almost massless up or down quarks and massive strange quarks. In that case, R ~ m^/(4pi?) ~ ml/{4:fi). 
The poles of the anomalous propagator are located at po = — i? ± (e^ -I- A^)^/^ — isgn(po)- As usual, we close the 
integration contour in the lower half plane. Let us denote the solution of the gap equation in the case i? = by Aq. 
Then, if i? < Aq, the pole with the positive sign of the square root is always included in the integration contour and 
we have 



A = 



47r2 



der, 



A 



(108) 



This result is, up to a small correction in the density of states that we have neglected here, identical to the gap 
equation for degenerate fermions, so A « Aq. If, on the other hand, R > Aq there only is a pole in the lower half 
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plane if ip > y/K^ — A'^. Carrying out the po integration again leads to the gap equation H108() . but with the Ip 
integration restricted by the condition just mentioned. This cuts out the infrared singularity at ep = and one can 
easily verify that the gap equation does not have a non-trivial solution for weak coupling. We thus conclude that a 
necessary condition for pairing is that 

ml < ml{crit) ~ 4pi?A(ms=0). (109) 



B. BCS theory: CFL phase 



So far, we have only dealt with a simple pair condensate involving strange and non-strange quarks. In practice, we 
are interested in a somewhat more complicated situation. In particular, we want to consider the transition between 
the color-flavor locked phase for small ms and the two-flavor color superconductor in the limit of large mg. This 
analysis can be carried out along the same lines as the toy model discussed above. We now consider the following free 
action il il 



S = 



d'^P / ,t i\ f iPo- ep)Xi - 2RX, AudXud + AusXus \ f ^ \ ^m^ 



(2 

where -0 is now a 9 component color- flavor spinor. Xi, X^, Xy^^ £^nd X^g are color- flavor matrices 

Xl = S^I'Sab, Xs = ^"'3^^3^^3 

-^ud — e e3a&, -^us — e ' t2ab + f- Elafc, 

where a, (3 are color, and a, 6 flavor indices. A^d is the gap for {ud) condensation, and A^s is the gap for {us) — (ds) 
condensation. Color- flavor locking corresponds to the case A^d — A^s, and the two flavor superconductor corresponds 
to = 0, And ^ 0. 

Flavor symmetry breaking is again caused by i? ~ ml/{App). The Nambu-Gorkov matrix equ. (|11U|I can be 
diagonalized exactly. The eigenvalues and their degeneracies are 



Po±{el + Al,f\ d = 3 

Po-R±(el + Aiy^\ d^2 



Po + i?±(e1 + AL)'^', d=2 
PO ± {el + 2Rep + 2R^ + 2AL + i A^, ±\sY'^ , d = 1, 1 

where 

1/2 



(112) 



(8AL(AL + 4i?,') + (AL-4i?(ep + i?))') . (113) 



The result becomes easier to understand if we consider some simple limits. If we ignore flavor symmetry breaking, 
i? = 0, and set Aud — A^s we find 8 eigenvalues po ± (e^ -I- A^)^/^ and one eigenvalue with the gap 2A. These 
states fill out an octet and a singlet of the unbroken SU{3)f symmetry of the CFL phase. If, on the other hand, 
we set Aus = we find 4 eigenvalues po ± (e^ -I- A^)^/^ while the other 5 eigenvalues have vanishing gaps. This is 
the spectrum of the Nf = 2 phase. If flavor symmetry is broken, we find that the SU{3) octet splits into two SU{2) 
doublets, one triplet and singlet. 

We note that in the presence of flavor symmetry breaking the first three eigenvalues, which depend on A^,; only, are 
completely unaffected. For the next 4 eigenvalues, which only depend on Aus , the energy pq is effectively shifted by R. 
This is exactly as in the simple toy model discussed above. It implies that for R > A^s , when we close the integration 
contour in the complex po plane, we do not pick up this pole. The last two eigenvalues are more complicated. They 
depend on both A^d and A„s, and they explicitly contain the flavor symmetry breaking parameter R. Nevertheless, 
the structure of the eigenvalues is certainly suggestive of the idea that for R < A^^, we have Aus — Aud, and the gaps 
are almost independent of R, while at i? ~ A°^ there is a discontinuity and Aus goes to zero. 

This is borne out by a more detailed calculation. For this purpose, we add a flavor and color anti-symmetric short 
range interaction 

C = ^{S^^S^' - 5"'6''''){SacSbd - SadSbc) (C'^2V'f) ^2^^^) ■ (IW) 



The free energy of the system is the sum of the quasi-particle contribution equ. H112|) and the mean field potential 

2 

'j^d 



V = 1/K ■ (A^^ -|- 2A^j,). There are two coupled gap equations, which can be derived by varying the free energy with 
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FIG. 12: Hierarchy of effective field theories in the CFL phase. 



respect to the two parameters lS.ud and A^^. A typical numerical result is shown in Fig. ^2 We observe that the 
flavor symmetry breaking difference A„(i — A^^ is quite small all the way up to the critical strange quark mass. At 
the critical mass, there is a discontinuous transition to a phase where A„s vanishes exactly. The value of the critical 
mass is very close to the estimate TOs = 2-\/^A0^. 

A number of authors have improved on the treatment pr esented in this section, in particular by studying the 
consequences of imposing electric and color charge neutrality jl02l llO.'j Il04j . In the BCS framework one finds that 
once charge neutrality is imposed the number density of up, down and strange quarks in the CFL phase is exactly 
the same, even in the presence of flavor symm etry breaking. As a consequence, the CFL phase does not require the 
presence of electrons to be electrically neutral 'IDS'] . 

Another interesting question concerns the possibility of additional phases that interpolate between the CFL and 
N f — 2 (2SC) superfluids. In |l06l Il07| it was shown that the charge neutrality constraint may help to stabilize 
gapless ('breached') CFL phases for TOj > ms{crit). Also, if the mismatch between the strange and non-strange Fermi 
surfaces is too large for BCS pairing to occur, pairing may still take place with a spatial ly varying sup erfluid order 
parameter [o^ . This phase is known as the LOFF (Larkin-Ovchinikov-Fulde-Ferell) phase |l08l Il09llll0| . The LOFF 
phase is distinguished by an interesting crystal structure 

All of these studies are based on an ansatz for the structure of the CFL phase in the presence of flavor symmetry 
breaking. This is somewhat unsatisfactory. In weak coupling and in the limit <C ms{crit) we should be able to 
perform rigorous calculations. Also, in the BCS approximation we find that the quark densities in the CFL phase 
remain exactly equal even if the strange quark mass or the electron chemical potential are non-zero. However, the 
CFL phase contains almost massless flavored Goldstone bosons, so the response to any external perturbation that 
can couple to Goldstone modes should not vanish. 

In the following sectio n we will show how to study the effect of a non-zero strange mass using an effective field 
theory of the CFL phase |ll2j | . This theory determines both the ground state and the spectrum of excitations with 
energies below the gap in the CFL phase. Using the effective theory allows us to perform systematic calculations 
order by order in the quark mass. 



C. CFL chiral theory 



For excitation energies smaller than the gap the only relevant degrees of freedom are the Goldstone modes associated 
with the breaking of chiral symmetry and baryon number, see Fig. 1121 The interaction of the Goldstone modes is 
described by an effective lagrangian of the form (112i] 



-t- 



J TT 



^Tr [VoSVoSt - vldm.^^] 



BTi{MT.'^) + h.c. 
AiTr(A/S^)Tr(MEt) + A2Tr(Af sUf S^) + A3Tr(AfEt)Xr(M^I]) 



(115) 



h.c 



Here E = exp(z(/)°A"//7r) is the chiral field, /tt is the pion decay constant and M is a complex mass matrix. The 
chiral field and the mass matrix transform as E LEi?^ and M LMR'^ under chiral transformations (L, R) G 
SU{?>)l X SU{i)f(. We have suppressed the singlet fields associated with the breaking of the exact C/(l)y and 
approximate U{1)a symmetries. As with ordinary chiral perturbation theory, the structure of the effective lagrangian 
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is entirely determined by the symmetries. At low density the coefRcients /tt, B, Ai, . . . are non-perturbative quantities 
that have to extracted from experiment or measured on the lattice. At large density, on the other hand, the chiral 
coefficients can be calculated in perturbative QCD. 

Superficially, equ. I|115|l looks exactly like ordinary chiral perturbation theory. There are, however, some important 
differences. Lorentz invariance is broken and Goldstone modes move with the velocity < c. The chiral expansion 
has the structure 

/:-/.^A^(|)'(|) {^n^r- (116) 

Loop graphs are suppressed by powers of p/(47r/^). We shall see that the pion decay constant scales as f.^ ^ pp. As 
a result loops are suppressed by p/pf whereas higher order contact terms are suppressed by p/A. This means that 
in the CFL chiral theory pion loops with leading order vertices are parametrically small as compared to higher order 
contact terms, whereas in ordinary chiral perturbation theory the two are comparable in size. 

Further differences as compared to chiral perturbation theory in vacuum appear when the expansion in the quark 
mass is considered. The CFL phase has an approximate {Z-2)a symmetry under which M —M and E — s- S. This 
symmetry implies that the coefficients of mass terms that contain odd powers of M are small. The (Zi)a symmetry 
is explicitly broken by instantons. The coefficient B can be determined from a weak coupling instanton calculation 
and B ~ {Kqcd/pfT, see Sect. HXEl 

A priori it is also not clear what the expansion parameter in the chiral expansion is. There are several dimensionless 
ratios that can appear, {m/pp), (m/A) and [m/ Kqcd)- The BCS calculations discussed in the previous section suggest 
that the CFL phase undergoes a phase transition to a less symmetric phase when m? / {2pp) ~ A. This result indicates 
that the expansion parameter is M'^ /{ppA). We shall see that this is indeed the case. However, the coefficients Ai of 
the quadratic terms in M turn out to be anomalously small. In Sect. HXDl we will show that 

/ AP\ 

A,M^ ~ A^M^ ~ f^A^ V~/ ' '■^^ 

compared to the naive estimate AiM"^ ~ /^A^[M^/(pfA)]. 

The pion decay constant /tt and the coefficients Ai can be determined using matching techniques. Matching 
expresses the requirement that Green functions in the effective chiral theory and the underlying microscopic theory, 
QCD, agree. The pion decay constant is most easily determined by coupling SU{Nf)L^fi gauge fields Wl,b. to the left 
and right flavor currents. As usual, this amounts to replacing ordinary derivatives by covariant derivatives. The time 
component of the covariant derivative is given by VqS = SqS + iWl^ — iTiWr where we have suppressed the vector 
index of the gauge fields. In the CFL vacuum E = 1 the axial gauge field Wp — Wr acquires a mass by the Higgs 
mechanism. From equ. (|115|l we get 

^-^liWL-Wnf. (118) 

The coefficients B and Ai^2.3 can be determined by computing the shift in the vacuum energy due to non-zero quark 
masses in both the chiral theory and the microscopic theory. In the chiral theory we have 



AS 



BTt{M) + Ai (Tr(M))^ + A2Ti{M^) + A3Tt{M)Tt{M^) + h.c. . (119) 



We note that as long as we keep track of the difference between M and il/t different O(M^) mass terms produce 
distinct contributions to the vacuum energy. This means that the coefficients Ai can be reconstructed uniquely from 
the vacuum energy. 



D. High density effective theory 

In this section we shall determine the mass of the gauge field and the shift in the vacuum energy in the CFL phase 
of QCD at large baryon density. This is possible because asymptotic freedom guarantees that the effective coupling 
is weak. The QCD Lagrangian in the presence of a chemical potential is given by 

C = i: {lip + ii-fo) ^ - i^LM^R - ^RM^i^L - \g%G%, (120) 
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where ^ + igA^ is the covariant derivative, M is the mass matrix and fi is the baryon chemical potential. If 
the baryon density is very large perturbative QCD calculations can be further simplified. The main observation is 
that the relevant degrees of freedom are particle and hole excitations in the vicinity of the Fermi surface. We shall 
describe these excitations in terms of the field V'+ {v, x) , where v is the Fermi velocity. The field V'+ (v, x) is defined on 
patches that cover the Fermi surface, see Fig. ^1 Soft collinear scatterings take place within a given patch whereas 
hard interactions can scatter Fermions from one patch to another. 

At tree level, the quark field ip can be decomposed as ip = ip+ + tp- where tp± = ^{l±d ■ v)'i/j. To leading order in 
1/pf we can eliminate the field using its equation of motion. For we find 

^-.L = ^ {id± ■ D^+,L + 70^^-+,/?.) • (121) 

There is a similar equation for ip-^R. The longitudinal and transverse components of 7^ are defined by ( 70, 7) || = 
iloi'vi'i-v)) and (7^)j. = 7^ — (7p)||. To leading order in 1/pf the lagrangian for the field is given by |ll3llll4lll5) 

- ((^i)' + MM^) + {R^L, M ^ Aft) + . . . ^ (122) 

with = (1, v) and i, j, . . . and a, 6, . . . denote flavor and color indices. In order to perform perturbative calculations 
in the superconducting phase we have added a tree level gap term 4''^ nCAai^bji'L r- the CFL phase this term has 
the structure Aai,bj — ^.{SaiSbj — SajSbi)- The magnitude of the gap A is determined order by order in perturbation 
theory from the requirement that the thermodynamic potential is stationary with respect to A. With the gap term 
included the perturbative expansion is well defined. 

The screening mass of the flavor gauge fields Wl^r can be determined by computi ng the co rresponding polarization 
function in the limit qo = 0, g — > 0. The relevant diagrams are shown in Fig. ^1 see |97llll7j| . The first two diagrams 
do not involve mixing between left and right-handed currents. The third diagram involves mixing between left and 
right handed currents and is unique to the CFL phase. We find 

Matching equ. against equ. l(TTH|l we get pTTllTTllTTgt 

21-81og(2) ( pI 



Repeating the matching calculation for the spatial components of the polarization tensor we get vl = 1/3 pTTj . 
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FIG. 14: Diagrams in the high density effective theory that contribute to the chiral flavor current polarization function at 
leading order in g. The green wiggly lines are flavor currents, while the black wiggly line is a gluon propagator. Squares are 
anomalous fermion self energy insertions, and the triangle is a gluon self energy insertion. 



Our next task is to compute the mass dependence of the vacuum energy. To leading order in 1/pf there is only 
one operator in the high density effective theory 

/: = (i^l+MM^i;L+ + ^i+M^Mi;R+) . (125) 

This term arises from expanding the kinetic energy of a massive fermion around p — pp. We note that MM^ /{2pp) 
and M'^M/{2pp) act as effective chemical potentials for left and right-handed fermions, respectively. Indeed, to 
leading order in the \/pf expansion, the Lagrangian equ. p22|l is invariant under a time dependent flavor symmetry 
ipL L{t)ipL, ipR R{t)^R where Xl = MMy{2pF) and Xr = M'^M/{2pf) transform as left and right-handed 
flavor gauge fields. If we impose this approximate gauge symmetry on the CFL chiral theory we have to include the 
effective chemical potentials Xl^r in the covariant derivative of the chiral field [97l |. 



(126) 



Xl and Xji contribute to the vacuum energy at 0{M'^) 

A£ = ^Tr [(MM^)(MUf) - {MM'^f] . (127) 

This result can also be derived directly in the microscopic theory '97]. The corresponding diagrams are exactly the 
same diagrams that appear in the calculation of /^r , Fig. ^1 but with the external flavor gauge fields replaced by 
insertions of equ. H125|l . We also note that equation (|127|l has the expected scaling behavior £ ^ /^A^[M^/(pi?A)]^. 

O(M^) terms in the vacuum energy are generated by terms in the high density effective theory that are higher 
order in the 1/pf expansion. These terms can be determin ed by computing chirality violating quark-quark scattering 
amplitudes for fermions in the vicinity of the Fermi surface [llg. Feynman diagrams for ql+Ql ^ Qr + Qr are shown 
in Fig. 115b . To leading order in the 1/pf expansion the chirality violating scattering amplitudes are independent of 
the scattering angle and can be represented as local four-fermion operators 

J0 = £r ((^f 'c^f Vv^gc^?)r-^^'"^ + (^f Vf )(^gVg)f^^^^) . (128) 

There are two additional terms with {L ^ R) and (M 4-> Aft). We have introduced the CFL eigenstates ifi"^ defined 
by = V^(A^)ai/\/2, A = 0, . . . , 8. The tensor F is defined by 

F^^^^ = i{Tr[A^Af(A^)^A^Af(A^)^] 

- ^Tr [A^Af(A^)^] Tr [A^Af (A^")^] |. (129) 
3 J 

The second tensor F involves both M and and only contributes to terms of the form Tr[Af Aft] in the vacuum 
energy. These terms do not contain the chiral field S and therefore do not contribute to the masses of Goldstone 



30 




FIG. 15: The left panel shows the diagram in the high density effective theory that contributes to the vacuum energy at order 
O(M^). The right panel shows how the effective four-fermion vertex can be matched against a chirality violating scattering 
amplitude involving a hard gluon exchange. 



modes. We can now compute the shift in the vacuum energy due to the effective vertex equ. H128|l . The leading 
contribution comes from the two-loop diagram shown in Fig. 1151 This diagram is proportional to the square of the 
superfluid density. We find 



3A2 

'4^ 



Tr[Af]) -Tr 



M ^ Aft 



(130) 



Note that a factor g'^ in the vertex equ. (|128|l cancels against a factor 1/g that comes from a logarithmic log(A) 
divergen ce in the superfluid density. Matching equ. H130|) against equ. H119|l we can determine the coefficients Ai.2,3. 
We find [nlllil 



3A2 
4^ 



^3 = 0. 



(131) 



We note that £ ~ /^A^(A/pi?) [M"^ /{ppA)] which shows that the coefficients Ai are suppressed by (A/pp). The 
effective lagrangian equs. (|122I128II can also be used to compute higher order terms in M. The dominant 0{AI^) 
term is the effective chemical potential term equ. (|127f) . Other 0{M^) terms are suppressed by additional powers of 



E. Instanton effects 



In the CFL phase spontaneous chiral symmetry breaking is dominated by order parameters of the form {{iptlj)^) 
and as a consequence the coefficient of the linear mass term in the chiral lagrangian is suppressed as compared to the 
quadratic mass term. Indeed, the CFL phase has an approximate (Z-2)a symmetry that forbids the linear mass term. 
The (Z-2)a symmetry is broken by instantons. This means that in the CFL phase the quark-anti-quark condensate is 
induced by instantons. 

We shall determine the coefficient B of the linear mass term by computing the instanton induced shift in the vacuum 



energy. The corresponding diagram is sho wn in Fig. 1161 In QCD with three flavors, the instanton induced interaction 

'2iVc + 1 



Qiagram is sho v^ 
between quarks is given by |l20l Il2ll Il22j| 

(27rp)6p' 



n{p,fi)dp 



2iVc 



8{Nc + 2) 



^L.,g2){'^R,h(T^,^i}L,g^) + {L ^ R) 



(132) 



Here, p is the instanton size, /i is the quark chemical potential, ,fi,gi are flavor indices and ct^^^ = ^[7^,71/]. The 
instanton size distribution n(p, fj) is given by 



n(p, 



C 



N 



87r2 
9^ 



2Na 



exp 



Stt' 

'gipy 



2 1 



exp [-A^/pV'] , 



(133) 
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FIG. 16: Instanton contribution to the 0{M) term in the vacuum energy. The diagram on the left shows the instanton term 
in the CFL phase at high baryon density. The squares are anomalous self energy insertions. The diagram on the right is the 
corresponding term in QCD at low baryon density. The squares are normal self energy insertions. 



Cat = 



0.466 exp(-1.679iVc)l-34^^ 



{N,-iy.{N, 

-61og(pA), 



2)! 

11 2 



(134) 
(135) 



At zero density, the p integral in equ. H133|) is divergent at large p. This is the well-known infrared problem of the 
semi-classical approximation in QCD. At large chemical potential, however, large instantons are suppressed and the 
typical instanton size is p ^ /^^^ <C A^^. We also note that at zero density the effective lagrangian equ. H132(l is 
derived by computing U{1)a violating Greens functions in the limit p ^ p~^. From a similar study at finite density 
one can show that a,t p ^ equ. (|132|l has to be int erpreted as an effective lagrangian for momenta in the vicinity of 
the Fermi surface |p| -pF < P"^ ~ M |80lll2l[l2^. 

To linear order in the quark mass one of the three zero modes is lifted. We find 



2(2^p)V^ 



n{p, p)dp ^^jy2 _ -^^ ^/i/2/3ggig2g3^-%33 



2iVc- 1 
27V, 



(V'fl,/iV'L,ffi)(V'i?.,/2'/'L,92) 
^ M\L^ R) 



(136) 



We can now compute the expectation value of equ. H136() in the CFL ground state _[6M l93l Il25| . To leading order in 
perturbative QCD we can use the mean field approximation {{ipip){'ip'ip)) ~ (ipip) (iptp) . The instanton contribution to 
the vacuum energy density is 



£ = 



16, 



4 „3 



'n'{.p,p)dp—{TTp) p 



9 



27r2 



Tr [M + M^] , 



where we have used the perturbative result for the diquark condensate in the CFL phase 

3V27r 



(137) 



(138) 



9 



V27r2 ; 



We note that for M = diag(m„, md, rus) the instanton contribution to the vacuum energy is indeed negative. Since the 
effective interaction involves both left and right-handed fermions the relative phase between the left and right-handed 
condensate in equ. (|138|l is important. Instantons favor the state with (V'lV'l) = —{^'R'i'R) which is the parity even 
ground state. Equation (|137|l for the vacuum energy can be matched against the effective lagrangian equ. (|119|) . We 
find 



B^C 



N- 



?7r4r(6) 
"3 ¥~ 
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27r2 
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9^ 



12 
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FIG. 17: Goldstone boson masses in the CFL phase as a function of the strange quark mass, from |97 



where we have performed the integral over the instanton size p using the one-loop beta function. The coefficient 
B is related to the quark-anti-quark condensate, (V'V^) = —25. We note that B is indeed parametrically small, 
B ^ {K/ppf. 

In QCD with three flavors chiral symmetry is broken both at small and at large density. The instantons contribution 
to the linear Goldstone boson mass term in the low density phase is shown in Fig. 116b . Again, the coefficient of the 
Tr(MS) term in the effective lagrangian is the instanton contribution to the quark conde nsate . There are strong 
arguments that instantons dominate chiral symmetry breaking in QCD at zero density, see |l22j for a review. This 
raises the question whether the instanton mechanisms of Fig. I16h and b are continuously connected. In a simple 
mean field calculation it was found that there is a first order phase transition that separates instanton induced chiral 
symmetry breaking at low and high density '93| , but this question certainly deserves further study. 



F. Kaon condensation 



Using the results discussed in the previous sections we can compute the masses of Goldstone bosons in the CFL 
phase. In Sect. HXCl we argued that the expansion parameter in the chiral expansion of the Goldstone boson masses 
is (5 = m? /{pf^)- The first term in this expansion comes from the O(M^) term in equ. (|115|l . but the coefhcients A 
contain the additional small parameter e = (A/pi?). In a combined expansion in 5 and e the 0{t5) mass term and the 
0{6'^) chemical potential term equ. p26|l appear at the same order. Instanton effects are suppressed by extra powers 
of [K/pp). To order 0(ei5, 6'^) the masses of the flavored Goldstone bosons are 



rriK. 



2pt 



2pf 



2pi 



— (m„ -I- mdjms 



-f2 

4A 

72 

./ TT 

J IT 



-,1/2 



1/2 



—mdirriu 



2 mu{md + iris) 
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Further studies of Goldstone boson properties in the CFL phase can be found in |ll5l Il26l Il27l Il28l Il29| . We observe 
that the pion masses are not strongly affected by the effective chemical potential fis = m1/{2pF) but the masses of 
the and are substantially lowered while the K~ and are pushed up. As a result the and meson 
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FIG. 18: This figure shows the phase structure of CFL matter as a function of the strange quark mass rris and the lepton 
chemical potential fiQ, from [q^I . 



become massless if 

m,|_, = 3.03.mf (141) 

see Fig. El For larger values of the kaon modes are unstable, signaling the formation of a kaon condensate. 

Kaon condensation implies that the CFL ground state is reorganized. For simplicity, we consider the case of exact 
isospin symmetry m„ — rud = m. Kaon condensation can be studied using an ansatz of the form S = exp(iQ;A4). The 
vacuum energy is 

Via) = -f^ i^!±-^^\iniar + (m°,)2(cos(a) - l)j , (142) 

where (m'j^)'^ ~ {4A/ f^)mu,d{mu.d + 'fns) is the 0{A'P) kaon mass in the limit of exact isospin symmetry. Minimizing 
the vacuum energy we obtain a = if m'^/{2pp) < to^ and cos(a) = (™/c-)^/Ms with — m1/{2pp) if /i^ > jti^^. 
In the kaon condensed phase SU{2)i x U{1)y is spontaneously broken to U{1)q. This coincides with the symmetry 
breaking pattern in the electroweak sector of the standard model. Kaon condensation also provides an interesting 
realization of Goldstone's theorem. Even though the number of broken generators is three, the number of Goldstone 
bosons is only two. This is related to the fact that one of the Goldstone modes has a quadratic dispersion relation 

The hypercharge density in the kaon condensed phase is given by 

ny (l - ^) . (143) 

This result is exactly analogous to the behavior of the charge density in the case of the dilute Bose gas, equ. H27|) . and 
Nc — 2 QCD, equ. (|101|l . We observe that within the range of validity of the effective theory, /is < A, the hypercharge 
density satisfies ny < Ap|,/(27r^). This means that the number of condensed kaons is bounded by the number of 
particles contained within a strip of width A around the Fermi surface. The maximum hypercharge density allowed 
by the effective theory corresponds to the case that essentially all strange quarks have been removed from the CFL 
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FIG. 19: Conjectured phase diagram of three flavor QCD with reahstic quark masses. 



wave function. This raises the question whether the CFL-2SC transition might be more comphcated than what is 
suggested by the simple BCS model discussed in Sect. 11X51 see Fig. ITTI 

In the limit of exact isospin symmetry the K° and condensed phases are degenerate. If charge neutrality is 
imposed then the K° condensed phase is preferred. The effective lagrangian equ. (|115|l can be used to study a number 
of question related to the structure and the low energy properties of the CFL phase. The phase structure as a function 
of the strange quark mass and non-zero lepton chemical potentials was studied by Kaplan and Reddy [o^, see Fig. 1181 
Other investigations have focused on the role of defect s such as v ortices (132.. il33j , as well as transport properties 
such as neutrino emissivity and thermal conductivity |l34i Il35i I13(t| . 



X. CONCLUSION: THE MANY PHASES OF QCD 



There are many issues in QCD at finite density that remain to be addressed. First of all it would clearly be of 
great help if some of the ideas discussed in these lectures could be tested using numerical simulations on the lattice. 
While some progress in this direction has been made, the most interesting regime of the phase diagram remains 
unexplored. Even more importantly, we have to find experimental or observational constraints on the properties of 
quark matter. In Sect. IVTl we mentioned the possibility of observing the QCD tri-critical in relativistic heavy ion 
collisions. Observational constraints on the properties of dense cold matter can be obtained from neutron stars. We 
refer the reader to for a recent review and a guide to the literature. 

On the more theoretical side we would like to improve our understanding of the phase diagram of QCD for realistic 
values of the strange quark mass. In Figs. I^IVl9l we have shown proposed phase diagrams for idealized versions of 
QCD. A conjectured phase diagram for QCD with three flavors and realistic quark masses is shown in Fig. ^] The 
reader will recognize many of the phases that we have discussed in the previous sections. However, from our discussion 
it is also clear that many of the phases, the transition lines, and the order of the transition are just guesses, based 
on mean-field or weak coupling arguments. Eventually, this problem will have to be resolved using observations and 
numerical simulations. In the mean time, however, we would also like to improve and systematize the theoretical 
approaches. This should be possible in both the limit of high and low density, using, respectively, weak coupling QCD 
and nuclear effective field theory. 
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